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Q-i Abstract: We consider the issue of vacuum stability and triviality bound of the 

singlet extension of the Standard Model (SM) with a singlet fermion dark matter 
(DM). In this model, the singlet scalar plays the role of a messenger between the 
SM sector and the dark matter sector. This model has two Higgs-like scalar bosons, 
and is consistent with all the data on electroweak precision tests, thermal relic den- 
sity of DM and its direct detection constraints. We show that this model is stable 
without hitting Landau pole up to Planck scale for 125 GeV Higgs boson. We also 
perform a comprehensive study of vacuum structure, and point out that a region 
where electroweak vacuum is the global minimum is highly limited. In this model, 
both Higgs-like scalar bosons have reduced couplings to the SM weak gauge bosons 
and the SM fermions, because of the mixing between the SM Higgs boson and the 
singlet scalar. There is also a possibility of their invisible decay(s) into a pair of 
DM's. Therefore this model would be disfavored if the future data on the (a ■ B)yy 
or (a ■ B)fj with V = 7, W, Z and f — b,r turn out larger than the SM predictions. 
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1 Introduction 

The Standard Model (SM) has been extremely successful in describing interactions 
between quarks and leptons down to ~ 10" 19 m, or up to a few TeV depending on the 
channels, and it would be complete if the Higgs sector is established experimentally, 
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clarifying the origin of electroweak symmetry breaking (EWSB) and masses for the 
SM chiral fermions and EW gauge bosons. 

However, the SM should be extended in order to accommodate 3 different direc- 
tions: 

• Neutrino masses and mixings cannot be explained within the context of renor- 
malizable SM. 

• Baryon number asymmetry of the universe requires a new source of CP viola- 
tion, beyond the CKM phase in the SM [1]. 

• Nonbaryonic cold dark matter of the universe should be included in the SM. 

Admittedly, the simplest and the most economic solution to the 1st and the 2nd 
problems is to invoke the seesaw mechanism [2] and leptogenesis by introducing extra 
singlet right-handed (RH) neutrinos [3]. This is also nicely fit to the idea of grand 
unified theory (GUT) based on SU(5) and SO (10). 

For the nonbaryonic cold dark matter, there are many candidates in particle 
physics models: axion and axino, the lightest supersymmetric particle (LSP) in SUSY 
models (neutralino or gravitino), the lightest Kaluza-Klein particle (LKP) in extra 
dimensional scenarios, to name only a few. Some of them are related with other 
problems in particle physics, such as strong CP problem or fine tuning problem of 
(Higgs mass) 2 , but there are many other models which are not related with other 
problems in particle physics. 

Another possibility is to rely on the principle of Occam's razor, namely the 
simplest extension of the SM with dark matter candidates. In terms of the least 
number of new degrees of freedom, a scalar DM model with Z 2 symmetry would 
be the simplest one. However the origin of Z2 symmetry is not clear, since it is 
put in by hand. The simplest DM model without ad hoc Z2 symmetry would be 
a singlet Dirac fermion CDM with conserved charge associated with a global dark 
U(l) symmetry. In Ref. [4], three of the present authors proposed such a scenario, 
by considering a Dirac fermion DM (-0) that couples to a real singlet scalar (S) (see 
also [5] and [6] for similar discussions). Writing the most general renormalizable 
lagrangian among these new fields (ip and S) and the SM fields, including the so- 
called Higgs portal terms, we could describe the DM physics (thermal relic density 
and direct detection), and collider phenomenology and electroweak precision tests 
(EWPT). Adding the singlet scalar improves the overall EWPT fits [4]. 

In this model, the Higgs phenomenology is modified in an important way by two 
different reasons: 

• There are two neutral Higgs-like scalar bosons, Hi and H 2 , which are two 
mixtures of the SM Higgs boson h and a singlet scalar s, with a mixing angle 
a. Couplings of Hi and H2 to the SM particles are reduced by cos a or sin a. 
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Therefore the production cross sections for if i=li2 's at colliders will be reduced 
by cos 2 a or sin 2 a compared with that of the SM Higgs boson with the same 
mass. 

• Both H\ and if 2 can decay invisibly into a pair of DM if kinematically allowed 
: Hi — > ipip. This would make more difficult to observe the H^s produced at 
colliders. 

These two independent mechanisms will make two Higgs-like scalar bosons have 
reduced signal strength a ■ B into specific final states [ see Eq. (3.3) ]. 

Recently, ATLAS and CMS reported a tantalizing hint for a Higgs-like boson with 
mass around 125 GeV [7-9]. Its couplings to the WW, ZZ and 77 are consistent 
with the SM predictions, albeit there are still large uncertainties because of limited 
statistics. More data accumulation is planned toward the end of this year, and we 
would learn much more about the detailed properties of the observed new Higgs-like 
boson. 

If the SM Higgs boson has mass around 125 GeV, the electroweak (EW) vacuum 
might be meta-stable or even unstable due to the quantum corrections from top quark 
loop [10-14] though large uncertainties in determining SM quantities including top 
pole mass do not allow to draw a firm conclusion on this issue [13, 15, 16]. Meta- 
stability might be still allowed as long as the tunneling time to wrong vacuum is 
longer than the age of our universe. However, if the primordial inflation is supposed 
to take place along the direction of Higgs field, for example as the case of Higgs 
inflation [17], the possible met a- stability of EW vacuum should be improved [14] 1 . 
Some new physics should be introduced well below the Planck scale in order to save 
this situation. It is the purpose of this paper to address this issue within the model 
proposed in Ref. [4]. In this model, there are only two more fields beyond the SM 
ones, the fermion DM ip and a real singlet scalar messenger S. Since 5* couples to the 
SM Higgs field directly, one can imagine that the EW vacuum in our model could 
be stable even if the new physics scale A is as large as Planck scale. In other words, 
the Planck chimney could be possible for Higgs mass around 125 GeV. 

However, once we introduce S and ijj into the SM, the vacuum structure can 
change significantly, which may give rise to various false vacua. For example, if S 
develops the vacuum expectation value (VEV), the Higgs potential in the S direction 
could take the form of a double-well potential. If it is tilted, it is no longer clear 
that the EW vacuum is the the global minimum. Furthermore, at the loop level 
if; contributes to the Higgs potential and may affect the vacuum structure as well. 
Nevertheless, such a vacuum analysis was often overlooked in the literature. 

In this paper, we investigate the vacuum structure and stability in the SM with 
S and ip. An effective potential approach is adopted to study the vacuum structure. 

1 Inflation due to Higgs field false vacuum [18-20] might be a possible alternative to Higgs 
inflation though the initial condition for inflation looks non-trivial to be realized. 
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We explore not only the EW vacuum but also other possible false vacua at the tree- 
and one-loop levels. At the tree level, we explicitly derive analytic expressions for 
the vacuum energies while the one-loop analysis exclusively relies on numerics. 

In order to examine the vacuum instability occurring at the high-energy scale, we 
use the renormalization group (RG) method. The /3 functions of all dimensionless 
couplings are derived at the one-loop level. As for the most relevant parameters 
such as the top quark Yukawa coupling, strong coupling and SU(2) doublet Higgs 
quartic coupling, we also include two-loop contributions coming from the SM sector. 
In addition to the vacuum stability, we also investigate the perturbativity of the 
quartic couplings up to the Planck scale. 

This paper is organized as follows. In Sec. 2, we describe the model and the 
relevant constraints from colliders and dark matter physics are discussed in Sec. 3. In 
Sec. 4, we discuss the vacuum structures, and various vacua are scrutinized carefully. 
The stability of EW vacuum taking account of the RG effects up to some new physics 
scale A is investigated in Sec. 5. The paper is summarized in Sec. 6. We collect 
the matching conditions used in our analysis in Appendix A, and the one-loop Higgs 
boson mass formulae as well as one-loop tadpole conditions are presented in Appendix 
B. 



2 The model 

We consider a SM gauge-singlet Dirac fermion DM (ip) with a real singlet scalar (S) 
that couples to the SM sector by the Higgs portal [4]. The dark sector is described 
by the lagrangian 



-•dark 



*l>{0 - m^W - XSijjip . (2.1) 



The most general renormalizable scalar potential including the Higgs portal interac- 
tions is given by 

V = -fi 2 H H^H + X H (tfH) 2 (2.2) 
+fi HS SH^H + l -\ HS S 2 H^H (2.3) 

+/4S + \™%& + l»'sS 3 + \\sS\ (2.4) 

where H is the SM Higgs field. 

In general, the neutral scalar fields develop nontrivial vacuum expectation values 
(VEVs), vh and vg- And we expand the neutral component of H and S as 

H=\^ {v l + h) \, S = v s + s, (2.5) 
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in the unitary gauge. Then the minimization conditions of the Higgs potential at 
VEVs give 

A4 = >^hv 2 h + H-hsVs + 2 X hsV 2 s , 

m s = (*s v s ~ A s % 5 ~ x hsv h . (2.6) 

We introduce the Higgs mixing angle a and the mass eigenvalues m i(=1 2 ) (mi < m 2 ) 
which diagonalize the Higgs mass squared matrix such that 



M 2 = | m L \ _ ( cos « sin « ] ( m i ) ( cosa ~ sin « j / 2 7) 
V Wj S J \ — sin a cos a J \ wij / \ sin a cos a 



The quartic couplings, A#, Xhs, Ass, in the Higgs potential can be expressed in terms 
of the Higgs mass parameters 



A 



n 



2 

lih 



2< 



*HS — • 



v s v H 

m ls + Vll V S + VHS} ,2 H / (2vg) - fj! s V S 

2vl 



X S = 7^2 ' ( 2 - 



S 

so that they are obtained as a function of mi,m2 and a which we take as input 
parameters. The mass eigenstates Hi (i = 1, 2) with masses are written in terms 
of the SM Higgs scalar h and the singlet scalar s as 

H i = h cos a — s sin a, 

H 2 = hsina + s cos a. (2.9) 

In the Higgs and dark sector we have 10 free parameters which are to be measured 
in the experiments: 

mi, m 2 , a, v H , v s , Us, a4> Hhs, wty(= m^ + Xv s ), A. (2.10) 

The constraints on these parameters from perturbative unitarity of electroweak gauge 
boson scattering amplitudes, EWPT, collider searches for Higgs boson(s), DM relic 
density, DM direct detection experiments, were given in Ref. [4]. And those from the 
vacuum stability and triviality will be considered in this paper. 

3 The constraints considered in Ref. [4] 

In Ref. [4] we considered the following observables which can constrain our model: 



- 5 - 



• the perturbative unitarity condition on the Higgs sector [21, 22], 

• the LEP bound on the SM Higgs boson mass [23], 

• the oblique parameters S, T and U obtained from the EWPT [24, 25], 

• the observed CDM relic density, Vt CI)M h 2 = 0.1123 ± 0.0035 [26], which we 
assume is saturated by the thermal relic ip, 

• the upper bound on the DM-proton scattering cross section obtained by the 
XENON100 experiment [27]. 

The first three conditions do not constrain the dark matter sector, and they are 
also relevant to the singlet scalar extension of the SM without dark matter. 

3.1 Perturbative unitarity of gauge boson scattering amplitudes 

The perturbative unitarity of scattering amplitudes for longitudinal weak gauge 
bosons in our model requires [21, 22], 

(m 2 ) = m\ cos 2 a + m 2 sin 2 a < ~ ( 70 o GeV) 2 , (3.1) 

3Gr p 

If mi ^ m2, Eq. (3.1) can be re-expressed as 

sin a 2 < - m\\ / (m 2 - mj) (3.2) 

which provides an upper-bound of the mixing angle as a function of ?772 for a given 
mi(« 125GeV). 

3.2 Collider bound 

In Ref. [4] we defined two ratios (i = 1,2) (what we called the reduction factor): 



~= a Ju B H su XsM (< = M), (3-3) 



where Xsu is a specific SM final state, which measure the reduced signal strength 
with respect to the SM. In terms of the SM Higgs ( singlet Higgs ) decay width p s ^ I ( hld ) 
with mass m^, i.e. without the effect of the mixing, we get 



„4 F SM 
_ c a L Hi 

Tl ~ „2FSM I „2r hid ' 

c a L Hi b a L Ht 

4 -pSM 



r 2 e 2 rSM , „2 Fhid j_ p ' (3-4) 

o; H 2 + C a L H 2 + 1 H % ^HiH x 
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We can see that the "reduction" of signal strength is a generic feature of this model, 
i.e. Ti < 1. If the future LHC data on n for some Xsm = VV,ff (V = 7, W, Z, 
f = b,r) is larger than 1, our model would be ruled out. In Ref. [4], we found that 
if r i ^ 0.7, we get r 2 < 0.2 for (m 2 > m^K 125GeV)). So the heavy scalar boson 
will easily evade the detection at LHC. 

If the Higgs splitting mode H 2 — > H\H.\ opens kinematically, it would provide a 
smoking gun signal for our model, like the four 6-jets, two-photons plus two 6-jets, 
four tau leptons, etc. The high-luminosity LHC machine can target these signals. 

3.3 The oblique parameters: S,T,U 

In our model the new scalar particle S can contribute to the W and Z boson self- 
energy diagrams, liww^zz, thereby changing the EWPT S,T,U parameters [4]. 
Explicit expressions for the oblique parameters in our model can be found in Ref. [4]. 
The result is that including singlet scalar improves the overall fit to the EWPT, 
which can be seen in Fig. 1. For m\ ~ 125 GeV, the mixing angle is constrained to 
be a < 0.4 when m 2 > 400 GeV. 

3.4 Dark matter relic density 

The observed DM relic density, ^cdm^ 2 ^ 0.1123 ± 0.0035 [26], is related with the 
thermally averaged annihilation cross section times relative velocity at freeze-out 
temperature roughly by 2 

ncDM „ „ 3X10-W 

The annihilation cross section of a DM pair is proportional to sin 2 2a. Since the 
EWPT and LHC observation of the SM-like Higgs boson restricts a to be small, the 
cross section is generically much smaller than is needed to explained the current relic 
density. This can be seen in Fig. 2 except for resonance regions. 

At resonance, if 7, = miYi/(Am^) 1, the non-relativistic approximation of the 
cross-section is [28] 



4vr 3/9 1/2 - X e B^l-B^il + Eif 12 Q 

KcsWlab)NR = — X J/ -7T i/2 7 ie Se « ^ ' K 6^) (3.6) 



where x = m^/T {z and = — 1 + m 2 /(4m 2 ). For example, ignoring H 2 — > H\H\ 
decay, we find 

pSMphid 
_i 
phi 



T.Bi (1 - B t ) = -pgM _i_ phid ' ( 3 - 7 ) 



2 There is a typo in this expression in Ref. [4] and we correct it here 
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Figure 1. The prediction of (S,T) parameters. We fixed the reference 

Higgs mass to be 120 GeV. The ellipses are (68, 90, 95) % CL contours from 

the global fit. The thick black curve shows the SM prediction with the Higgs 

boson mass in the region (100,720) GeV. The red, green dots correspond to 

a = 45°, 20°, respectively. The dots are for the choices (mi, m,2)(GeV) = 
(25, 125), (50, 125), (75, 125), (100, 125), (125, 125), (125, 250), (125, 500), (125, 750) from 
above for each color. 



If m 2 decays dominantly to dark matter, 



3y / 2 

T H2 B H2 (1 - B H2 ) ~ sin 2 a G F m 2 f m 2 (3.8) 

3 2 f m t\ 2 

= — sin a — m 2 , (3.9) 

87T \ V / 



and hence 

2 



(a resUlab ) NR ~ ^^e-^ sin 2 a ( (3.10) 



6 x 10- 5 sin 2 a ( GeV" 2 (3.11) 

V m 4> J 

where we have used x = 25 and 6j = 1/x in the second line. Hence, even if a might 
be constrained to be small, a right amount of CDM relic density can be obtained as 
long as is in the band of s-channel resonance. 
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Figure 2. Dark matter thermal relic density (f^cDM^ 2 ) as a function of mi for m\ = 
125 GeV, A = 0.4, a = 0.1 and m f = 100, 500, 1000, 1500 GeV from top to bottom at right 
side. The dotted red line corresponds to the observed value, ^cdm^ 2 — 0.112. 

We used the micrOMEGAs package [29] for numerical calculation of DM relic 
density and direct detection cross section. In Fig. 2, we show the CDM relic density 
as a function of m 2 for various choices of = 100, 500, 1000, 1500 GeV, with A = 0.4 
and a = 0.1. We can always find out the rri2 value which can accommodate thermal 
relic density of the singlet fermion CDM ip. Note that there is no strong constraint 
on the heavier Higgs with a small mixing angle a, because H 2 would be mostly a 
singlet scalar so that it is very difficult to produce it at colliders, and also it could 
decay into a pair of CDM's with a substantial branching ratio. 

3.5 Direct detection 

The null searches of the dark matter-proton scattering puts strong bounds on its 
spin-independent (SI) cross section [27]: 

a p < 2 - 10 (l(T 9 pb) (3.12) 

for the CDM in the mass range mcDM = O(10 — 100) GeV. The spin- dependent 
scattering cross section is zero in our model, because the scattering is due to two 
Higgs-like scalar bosons. Since the a p is proportional to the (cr ami v)f z , the large 



- 9 - 



annihilation cross section need for the relic density would also give large DM-proton 
scattering cross section, violating (3.12). 

As pointed out in [4], there would be a destructive interference between Hi and 
H 2 contributions to the scattering amplitude due to orthogonality of the Higgs mixing 
matrix, which is a very generic aspect in case there are extra singlet scalar bosons 
that can mix with the SM Higgs boson [30]. Hence, for regions m 2 — mi <C m 1; a 
cancellation occurs in a p and even the large A, a regions are only weakly constrained 
(see Fig. 4 of Ref. [4]). We also note that a p and (a arm v){ z are not strongly correlated 
near the Higgs resonance where the relic density can be explained. This helps to 
evade the strong bound on a v while accommodating the correct CDM density in the 
universe. This opens up a very interesting parameter space for Higgs boson search 
at the LHC, making one or two of the Higgs-like scalar bosons can decay into a pair 
of DM's with a substantial invisible branching ratio(s). 

3.6 Comparison with the effective lagrangian approach 

In this subsection, we would like to compare our model with the so-called Higgs 
portal fermion dark matter model [6, 31], where the singlet scalar 5* is presumed to 
be integrated out, resulting in the following model lagrangian: 



Within this model, there is only one Higgs boson and its coupling to the DM is 
strongly constrained by the direct detection experiments. This result is very differ- 
ent from our analysis [4], where there is a generic cancellation between Hi and H 2 
contributions in the direct detection rates. In fact, o"si depends also on (sin a cos a) 2 , 
and it becomes zero when we ignore the mixing between the SM Higgs boson and the 
singlet scalar S (see Eq. (3.16) of Ref. [4]). This result can never be obtained in the 
approach based on the above effective lagrangian (3.13). In our case the correlation 
between Hi — ip — ip coupling and the direct detection cross section is not that strong 
compared with the results in Ref. [31]. It is important to consider the renormaliz- 
able models in order to discuss phenomenology related with the singlet fermion dark 
matter and Higgs bosons. 

The same arguments also applies to the Higgs portal vector DM models, which 
is assumed to be described by the following lagrangian: 



Although this lagrangian looks power-counting renormalizable, it is not really renor- 
malizable. This is well known from the old intermediate vector boson theory for 
weak gauge boson W ± . In order to give a mass to a spin-1 gauge boson, we need 
some symmetry breaking agency. Assuming a new complex scalar (px breaks the 




(3.13) 



c = -m 2 v v^ - -^-rfHv,y» - -fiy^f . 



(3.14) 
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gauge symmetry spontaneously, one ends up with a new scalar boson from cf>x which 
would mix with the SM Higgs boson by Higgs portal. Therefore there will be two 
Higgs-like scalar boson in the end, and phenomenology in the scalar sector should 
be similar to that of the model described here and in Ref. [4]. We leave the detailed 
discussions of this issue for the future publication [30]. 

4 Vacuum structure 

Because of the presence of the singlet scalar, the vacuum structure of this model is 
not that trivial. Since the Higgs potential is the quartic function of the Higgs fields 
(at the tree level), there could be another nondegenerate local minimum in the singlet 
Higgs direction unless some symmetry exists. If that is the case, our EW vacuum 
may not be global and its stability is unclear. In addition to this, the EW vacuum 
could be destabilized at a high energy scale by the RG effect of the top quark as in 
the SM. We separately examine the vacuum stability at the EW scale and the high 
energy scale. In this section, we focus on the former, and the latter will be discussed 
in the next section. 

4.1 Tree level analysis 

Let us first consider the vacuum structures of our model at tree level 3 . In this 
analysis, \j? H and m| are determined by Eq. (2.6) with fixed vh and v$- The tree- 
level effective potential then takes the form 

rr I \ , 4 2 2\ , ^HS (2 2 ^- V H l Ps\ 

V {(p H , <Ps) = -£VPh- 2v h¥h) + -y- I <Ph<Ps - ¥h v s - 7, J 

, ^HS , 22 22 2 2\ , 3 f 1 

+ —{<Ph<Ps ~ Vh v s ~ v h 1 Ps) + yPs - 2~ J 

+ f (f a - ftwl) + ^V 5 - 2«M), (4-1) 

where <pn and tps are constant background fields. To avoid the potential unbounded 
from below, we impose 

X H >0, A s >0, X 2 hs <4X h Xs, (4.2) 

where the last condition is needed for Xhs < 0. 

Unlike the SM, there is a possibility that Vq((Ph, fs) has a global minimum which 
is different from the prescribed vacuum (vh,vs)- Following the Refs. [33, 34], we 

3 The vacuum analysis of the singlet extension of the SM within the electrowcak phase transition 
context can be found in [32]. 
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define the various vacua as follows: 



EW : v H = 246 GeV, v s = vf, 
SYM : v H = v s = 0, 
I-.vh — O, v s ^ 0, 
II : v H ^ 0, v s = 0, 
III : v H ^ 246 GeV, v s ^ v$, 



(4.3) 
(4.4) 
(4.5) 
(4.6) 
(4.7) 



where v 1 ^ is the prescribed v$- In the phase III, although both v h and vs are nonzero, 
they are different from the prescribed vacuum (EW phase). Those various vacua are 
found by solving the following equations 



dV 



dip H 
dV 



d<p< 



LP S =V S 



^HV 2 H + /J-HSVS 



A 



HS -2 2 

7T v s ~ /% 



0. 



(4.8) 



X s v s + n' s v 2 s + [ m 



«5 + ^«i + /4 = 0. (4-9) 



Note that one of the solutions corresponds to the EW phase. For the EW phase to 
be the global minimum, we require 



V (vh,v s ) < V Q (v H ^ v H ,v s ^ v s ) 



(4.10) 



where vh,s denote the VEVs in the SYM, I, II and III 4 . To begin with, we demon- 
strate a comparison between the EW and I phases. The vacuum energies of the both 
phases follows. 



V f 



(EW), 



A 



- 'H 4 

[VH: V S ) = —V H - — V H V S ~ 

T /(I)/ n - \ V'S- Vs-3 ^S_4 

K (0, vs) = —v s - —v s - —v s . 



Vhs 



^tW + f 8 "^-^ (4 ' li: 



/4, 



2 6 ° 4 

Here, we define A^ I_EW )Vo by taking the difference of the two vacuum energies 



(4.12) 



A P-*wv = y o m(o 



v,s) 



Vr 



(EW) 



A// 4 fiHS 2 
Vir H "' 

4 * 4 



A 



. "#S 2 2 i A^A' /- 
u tf u S + —1%V S + -j( V S ~ ^s) ~ 



VH,V S ) 

t4 



if' 





(4.13) 



To satisfy the condition (4.10), A^ I_EW ^Vo should be positive. However, it could be 
negative if vs gets large, which we will illustrate in the following. For simplicity we 
take us — 0. From Eq. (4.9), it follows that 



v s = 0, 



2A 



s 



(4.14) 



4 Since we will not consider a case in which both ^lhs an d A*s are simultaneously zero in the 
following discussion, the II phase would not be realized. 
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where the second solution corresponds to the I phase, and > A\ s m 2 s should be 
satisfied for real solutions. The vacuum energy of the I phase is reduced to 



m s + 



6Ac 



(4.15) 



Therefore, for the large values of fj/ s and m 2 s with their appropriate signs, we may 
have A( I-EW )Vo < 0. As we will see later, the large m 2 can induce such a case. 

Now we consider another solution in Eq. (4.8). The nonzero vh is expressed as 



V 2 H 



A 



H 



A 



HS -2 



S 



/4 



(4.16) 



where the real solution of vh enforces ii 2 H > Xhs^s/^ + Vhs^s- Plugging this into 
Eq. (4.9), we have a cubic equation for vs- If the cubic equation has only one real 
solution, it is nothing but v$ and the III phase cannot exist. On the other hand, if 
the cubic equation has other real solutions, and simultaneously Eq. (4.16) has a real 
solution, the III phase would appear. In such a case, it should be checked whether 
the energy level of the EW vacuum is lower than that of the III phase. Let us denote 
the real solutions other than v$ by and v^ 2 ' 
gives a local maximum and thus v 
III phase. 

Before going to the numerical analysis of the vacuum structures, we will obtain 
a range of m 2 consistent with the global vacuum conditions. To make the analysis 
simpler, we set a — — 0. In the limit of y/\ns v H <^ \/\s\ v s\, m 2, we may find 



— (2) 

s aii^ u s . We define v s as the solution that 
s and its corresponding solution of vh yield the 



-2X s v s + — , ml ~ X s v 2 s - ml. 



The vacuum energy of the EW phase then can be cast into the form 



(EW) / s 

o (v H ,v s ) 



A 



Ha 



V's-i A 



v H —v Q - 



A 



v s ~ 



■ 2 2 
:V s m 2 



A 



H 4 



(4.17) 



(4.18) 



Requiring A^ SYM - EW ^V = V (SYM) (0,0) - V^ EW) (v H ,v s ) > yields the lower bound 

of 777.2: 



A> 




vs\ < m 2 . 



(4.19) 



Furthermore, in the limit of v^-VsK'sl < m 2, A 1 - 1 EW ^Vo can be approximated as 



A (I-EWW 



777 



2\4 



12 A|4 



+ 77 



1 (m|) 3 



6 Af 7,1 



(4.20) 
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Figure 3. (Left) The VEVs of the prescribed vacuum (EW phase) and the nontrivial 
vacuum (III phase) as a function of a. (Right) The vacuum energies of the two vacua 
as a function of a, where V = Vq(vh,vs) an d Vq IU ^ = ^b(^#\ ^g^)- Here we take 
mi = 125 GeV, m 2 = 200 GeV, v s = 500 GeV, X HS = 0.01, A 5 = 0.2, fi S = 0. 



Therefore, A^ 1 EW )Vo > gives the upper bound of 777.2: 

m 2 < y/2X^\v s \. (4-21) 

Similarly, we may also obtain another constraint on m 2 from A( m - EW V > 0. In- 
stead of doing so, we investigate the vacuum structures numerically. 

We begin with a case in which the EW vacuum becomes the local minimum and 
the III phase can be the global minimum. The representative example is shown in 
Fig. 3. We here take m 1 = 125 GeV, m 2 = 200 GeV, v s = 500 GeV, \ HS = 0.01, A 5 = 
0.2. In the left panel, the red lines represent the prescribed VEVs (vjj,vs) =(246 
GeV, 500 GeV) and the upper and lower blue dotted curve denote and Vg 
which significantly depend on the values of a, and v s is mostly negative. In the 
right panel, the vacuum energies of the EW and III phases (V^ EW ^ andV^ 111 ^) are 
shown, where Vq 111 ^ is given by Vo(ujy , Vg ). We can see that Vq E < Vq holds 
only up to a ~ 0.3 rad, and beyond this, the III phase becomes the global minimum. 

As suggested by Eqs. (4.19) and (4.21), the region where the EW phase is the 
global minimum is also highly limited by the values of m<i for a given v$ and \g. In 
Fig. 4 we illustrate such constraints. In the upper panel, the vacuum structure is 
shown in the a-m-i plane, taking the same input parameters as those in Fig. 3. It is 
found that the EW phase can be the global minimum only for 160 GeV< m 2 < 320 
GeV, and in most of the parameter space, the III phase is the global minimum. It 
should be emphasized that the possible range of m 2 around a ~ rad is completely 
consistent with the analytic formulae (4.19) and (4.21). From this observation, al- 
though we have not worked out the analytic formula from A^ II_EW ^Vo > 0, the mass 
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Figure 4. The tree-level vacuum structures in the a-mi plane. (Upper) vs = 500 GeV; 
(Lower Left) v s = -500 GeV; (Lower Right) v s = -1000 GeV. 



bounds on the second Higgs boson is more or less the same as (4.19) and (4.21). 
In passing, we also find a region where the I phase becomes the global minimum at 
a ~ rad and > 350 GeV. 

If the sign of Vg is changed, namely, Vg = —500 GeV is taken but keeping the 
rest of the input parameters, the vacuum structure is drastically changed as shown 
in the left lower panel. In this case, depending on a and 7712, the EW, SYM, I and III 
phases can become the global minimum. From this plot, we obtain the upper bound 
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of the heavy Higgs boson mass, i.e., m 2 < 500 GeV around a = 0.8 rad. Since a is 
relatively large, this upper bound cannot be obtained from (4.21). 

In the right panel, we set vs = —1000 GeV. We can see that the global minimum 
region of the EW phases is significantly affected by the value of vs, rendering m 2 be 
as large as 800 GeV around a = 0.65 rad. The possible lower value of m 2 is also 
pushed upward, m 2 > 300 GeV. For the small a region, the allowed range of m 2 
nicely agrees with the analytic formulae (4.19) and (4.21). 

4.2 One-loop level analysis 

The tree-level Higgs potential receives the quantum corrections from the particles 
that couple to the Higgs fields. Therefore, the vacuum structure may change at the 
loop level. In our model, the DM in the hidden sector can also affect the Higgs 
potential through its couplings to the singlet Higgs field. So far, such an effect on 
the vacuum structure has not been investigated in the literature. 

We analyze the vacuum structure using the one-loop effective potential [35] 

^iWs) = E^(^f ( 4 - 22 ) 

which is regularized in the MS scheme; c = 3/2 for scalars and fermions and c = 5/6 
for gauge bosons, fi is a renormalization scale which will be set on m 2 . m is a 
field-dependent mass, and % denotes particle species which is explicitly given by 
% = H12, G°, G ± , W, Z, t, b, and their degrees of freedoms (n*) are respectively given 
by 

n Hl = n H , 2 = n G o = 1, n G ± = 2, n w = 6, n z = 3, 

n t = n b = -12, = -4. (4.23) 

Unlike the tree-level analysis, it is impossible to obtain the analytic formulae of the 
vacuum structures. We thus numerically minimize 

V eS ((p H , Vs) = V (tp H , (p s ) + V^ifn, tp s ) (4.24) 

and find a global minimum. Since we are considering the vacuum stability at the low 
energy scale, we will concentrate on the fields space below 10 TeV. 
We impose the one-loop tadpole conditions as 

/dV cS \ /dV \ ^ m? /dfri%\(, mj 1\ n , Ane . 

where = fn, Vs, an d (' ' ' ) is defined such that a field-dependent quantity is 
evaluated in the EW vacuum, namely, <pn = Vh and tps = vs- Since we are taking 
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vh, v$, mi, m 2 and a as the input parameters in place of fi 2 H , m|, X H , u>hs an d 
we solve the following coupled equations numerically 



vh \ dip H J v s \ dips J 



1 /dV eS \ 1 /8V eS \ 



(4.26) 



2 2 2 2-2 n 

2 2-2 2 2 n 

m ss — m 1 sm a — m 2 cos a = 0, 
m L + ( m i — m 2) s i n a cos a = o, 



(4.27) 
(4.28) 
(4.29) 



where m\ h , m 2 ss and m\ s are defined by the mass matrix of the Higgs bosons at 
the one-loop level. The explicit expressions of one-loop quantities are listed in Ap- 
pendix B. Since the treatment of the Nambu-Goldstone (NG) boson loop contribu- 
tions are somewhat tricky and numerically unimportant, we will not take them into 
account in the analysis here. 

In Fig. 5, the vacuum structures at the one- loop level are shown. The input 
parameters in the Higgs sector are the same as in Fig. 4. As for the parameters in 
the DM sector, as an example, we set 



In the upper plot, the significant difference between the tree and one- loop results is 
not observed in the entire region. As for the lower two plots, on the other hand, 
the region of the SYM phase appearing in the tree-level analysis vanishes, and the 
I phase region is enlarged instead. The reason is the following. At the tree level, 
(vh,vs) = (0,0) solution can exist if [is = as can be seen from Eqs. (4.8) and 
(4.9). At the one- loop level, on the contrary, there remain the constant terms in the 
tadpole condition for ip s even after taking v H = as we can observe in Eq. (B.7). 
Those constant terms are proportional to fins or m ^ which are always nonzero for 
the input parameters we are choosing here. Therefore, the SYM phase can never be 
realized at the one-loop level. 

For i>s = —500 GeV, the EW phase region is virtually unchanged and thus the 
allowed range of remains the same. For vs = —1000 GeV, on the other hand, 
although the vacuum structure in the small a region is not much affected by the 
one-loop contributions, the EW phase region is significantly distorted in 0.6 rad < 
a < 0.8 rad, and the maximal rri2 can reach 1000 GeV around a = 0.8 rad. 

Here, we remark that since the chosen can be the half of r«2 in the allowed 
region, the DM relic density can be explained by the resonance effect as demonstrated 
in the previous section. 

In order to see the DM loop effects on the vacuum structure, we vary A as 




A = 0.2. 



(4.30) 



A = ±0.01, ±0.6, 



(4.31) 
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Figure 5. The one- loop level vacuum structures in the a-m.2 plane. (Upper) vs = 500 GeV; 
(Lower Left) v s = -500 GeV; (Lower Right) v s = -1000 GeV. We take = y/X^\v s \/2. 



taking m^ = \vs\ and the rest of the parameters are the same as the previous cases. 
For the moment, we do not take account of the DM relic density constraint. Our 
findings are shown in Fig. 6. In the upper (lower) panels, v$ = —500 (—1000) GeV is 
taken, and in the left (right) panels, we set A = 0.01, 0.6, (—0.01, —0.6), where the red 
curves correspond to A = ±0.01 and the blue dotted curves represent A = ±0.6. For 
the positive A, the primary effect due to the increase of A is the shift of the EW phase 
region to the right side, correspondingly, the minimal values of gets enhanced by 
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Figure 6. The effects of A on the vacuum structures. The red straight curves corresponds 
to the case of |A| = 0.01 and the blue dotted curves denotes to the case of |A| = 0.6. (Upper) 
v s = -500 GeV; A = 0.01,0.6 (left) and A = -0.01,-0.6 (right). (Lower) v s = -1000 
GeV; A = 0.01,0.6 (left) and A = -0.01,-0.6 (right). 



about 25 (50) GeV around a ~ rad, and the maximal values of m 2 is enhanced by 
about 50 (200) GeV around a ~ 0.8 rad in the case of v s = -500 (-1000) GeV. 

For the negative A, on the other hand, the EW phase region is more sensitive to 
the change of A, especially vs = —500 GeV shown in the upper right panel. 

The III phase region which can exist in the A = —0.01 case is gone if A = —0.6, and 
the EW phase region is much enlarged. In this case, the range of mi < m 2 < 200 
GeV is allowed for any a G [0, 1]. The maximal value of m 2 , however, is reduced to 
400 GeV which is realized at a ~ 0.75 rad. 

In the case of vs = —1000 GeV, the minimal value of m 2 is increased by about 
100 GeV, and the maximal value of it gets bigger by about 150 GeV when A is 
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changed from —0.6 to —0.01. 

Here, we comment on the DM relic density constraint. The choice of m^, = \vg\ 
with Eq. (4.31) does not yield = 7712/2 in the EW phase region. Therefore, the 
right amount of the DM relic density may not be guaranteed. It turns out that the 
change of A has little effect on the vacuum structures once we take = 7772/2. 

In summary, in this section we showed that the diverse types of the (false) 
vacua are realized in this model, which is due to the presence of the singlet scalar 
field. Most important consequence is that the EW vacuum is not always the global 
minimum and often becomes the metastable state. As explicitly demonstrated here, 
the global minimum condition for the EW vacuum can eliminate the large portion 
of the parameter space. As a result, we can obtain the strong bounds on m 2 . For 
a < 0.2 rad, we find \/\s/2\vs\ < < V^s\ v s\- Is should be stressed that this 
mass bound exclusively depends on vs and As and not on the value of \hs- For 
v s ^ v Hi however, the above 7772 range would not be valid any more. The EW phase 
regions appear somewhere near a = ir/2 in which 777 2 becomes the SM-like Higgs 
boson. 

At the loop level, the DM also contributes to the effective potential. It is found 
that the vacuum structure has some sensitivity to the magnitude of A, rendering the 
viable ranges of a and 777 2 changed. However, such effects would be diminished once 
777^, = 777,2/2 is imposed, which is indeed the case if we wish to saturate the observed 
thermal relic density of CDM by a singlet fermion CDM. 

5 Vacuum stability 

Recent results of ATLAS and CMS experiments may indicate rrih ~ 125 GeV. If 
SM is assumed to be valid up to a very high energy scale, for example, GUT or 
Planck scale, the SM Higgs potential realizing such a light higgs faces the problem 
of vacuum instability, that is, the existence of non-SM deeper vacuum or unbound- 
from-below caused by the negativity of quartic self-coupling of Higgs at large field 
region. However, a simple extension of SM changes this situation drastically. As 
shown in the previous section, including a singlet which couples to SM via Higgs 
portal make vacuum structure very complicated. In this section, we show how the 
SM picture of instability problem of Higgs potential at high energy scale is changed 
in our model. 

Before getting into the analysis, we fix the low energy boundary quantities at 
the scale of Z-boson pole mass. 



a em (M z ) 



127.926 



1 



a 2 (M z ) 



a cm (M z ) 
sin 2 W 



a 3 {M z ) = 0.1184 



(5.1) 



where 



M z = 91.188 GeV 



sin 2 9 



•w = 



0.2312. 



(5.2) 
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We run up the couplings to the scale of top-quark running mass in MS-scheme. 
Denoted as m t , the running mass is obtained from the well-known formulas listed in 
Appendix. For the top-quark pole-mass [36] 



M t = 173.2 GeV, 

we find 

m t (m t ) ~ 164.0 GeV. 
Correspondingly, the running top- Yukawa coupling is given by 

X t {m t ) = V2m t (m t )/v. 



(5.3) 



(5.4) 



(5.5) 



Note that choosing M t as the matching scale results in m t (M t ) ~ 163.5 GeV which 
is smaller than m t (m t ) by about 0.5 GeV. As shown in Fig. 7, this difference results 
in about 0.5 GeV higher instability scale of SM Higgs potential. It is similar size to 
the uncertainties of the top-quark pole mass and associated instability scale. 

In order to obtain the running Higgs quartic coupling at m t (m t ), we solved 
Eqs. (4.26)- (4.29) numerically, ignoring momentum dependent corrections which is 
about %-level contributions [37] 5 . 

5.1 Tree-level (mixing) effect 

The first thing we should note in a extension of SM like ours is how the mass of the 
SM-like higgs is determined. For the fixed VEV of SM Higgs field, Xh is no more 
the only parameter which determines the mass of SM-like higgs. As studied in Refs. 
[38, 39], tree-level contribution of a singlet scalar can remove the instability problem 
in a very simple manner. Similarly, the tree-level effect of our model on the stability 
of SM Higgs potential can be read off from the mixing effect on the SM Higgs. At 
tree level, the SM Higgs quartic coupling can be written as 



A 



H 



1 + tan (a) — ^ 

m{ 



cos [a] 



mi 



(5.6) 



where a is the mixing angle. If a = 0(no mixing), one obtains Xh = A)^ 1 = m\j (2v 2 ) 
with mh = 125 GeV. However, once the mixing is turned on, Xh can be much larger 



than A| M . 



It may turn out that Higgs is very SM-like, so the mixing between the singlet and 
the SM Higgs might have to be highly suppressed. Even in such a case, we can still 
have a sizable increase of Xh to achieve the stability of Higgs potential by pushing 
up m2 to its unitary bound at most. If small, the mixing angle is given as 



a 



m hs 




vh (^hs + XhsVs) 


m ls - m lh 




m2 ss ~ m lh 



(5.7) 



5 In case of SM, a %-level correction to Higgs mass results in about an order of magnitude change 
of the instability scale. 
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Figure 7. The scale dependence of SM Higgs quartic coupling Xh for = 125 GeV. 
The solid and dashed red lines correspond to the cases with MS top- Yukawa couplings 
obtained at the matching scale mt(Mt) — 163.5 GeV and mt(mt) — 164 GeV, respectively. 
The instability scale is different from each other by about 0.5 GeV. 



In the limit of large v$ where the heavy singlet-like Higgs is likely to be decoupled, 
Xh is approximated as 

Xh - XT + ^ (5.8) 

which is the same as the case of Ref. [38, 39]. On the other hand, if hhs ~ ~Xhs v Si 
the tree-level effect on Xh is negligible, but the loop-effect from extra particles could 
be still large enough to remove the vacuum instability as long as Xhs is sizable. 

In Fig. 8, we show Xh contours in (m 2 ,a)-plane. If the RG-running of Xh is pure 
SM-like, vacuum stability requires m h > 130 GeV (corresponding to Xh > 0.139) at 
2-loop level. We can see in the figure that, even if a < 0.1, vacuum can be stable if 
m 2 > 400 GeV. 
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Figure 8. Tree-level mixing effect on A# for mi = 125 GeV in (m2,a) plane. The x-axis 
(a = 0) corresponds to pure SM where Aij(mj) = 0.1277. In SM, vacuum stability up to 
Planck scale requires Aff(mt) > 0.139(white region) corresponding to rrih ~ 130 GeV for 
Mt = 173.2 GeV and a s (Mz) = 0.1184. The dashed blackdine corresponding to a = 0.1 is 
a reference line that might be imposed from LHC Higgs searches. 



5.2 Loop effect 
5.2.1 RG equations 

If a is limited to be small due to constraints from collider experiments and/or direct 
searches of dark matter, we may have to resort to the additional loop contributions 
coming from the Yukawa couplings, Xhs, an d A. 

The /3-function of a coupling A, in the renormalization group equation is defined 

as 

fi x . = dX/dlafx (5.9) 
where [i is the renormalization scale. For SM gauge couplings, the /^-functions are 
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given as 



1 



9a 167r 2^" ' (-^2) 

where, for notational convenience, we have redefined the SM couplings as 

9i = 9 , 92 = 9 , 93 = 9s 



and 



41 19 
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(5.10) 
(5.11) 

(5.12) 
(5.13) 

(5.14) 



For dimensionless couplings (including the top Yukawa coupling) in the scalar po- 
tential, the 1-loop /3-functions £1X6 cLS follows. 
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For the top- Yukawa and quartic self couplings of Higgs field, the 2-loop /3-function 
contributions are as follows. 
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Note that /3\ HS is linearly proportional to Xhs, hence Xhs does not change its sign 
during its RG-running. A does not change its sign, too, and it increases or decreases 
monotonically, depending its sign. A appears in the RGEs of Xhs an d Xs as a squared 
one, so its sign does not affect RG-running of those couplings. 

As can be seen from Eq. (5.19), the RGE of A does not depend on any other 
couplings at 1-loop level. Hence it can be solved easily, giving the solution 



A(/i) 

We find that /3 A < 1 if 



which corresponds to 
X(m t ) 



In" 



A 2 (/x ) 8vr 2 ji Q 



-1/2 



(5.22) 



1/3 



A(M P1 ) < ( i^l ) (5.23) 



In 



A(Mpj) 2 8tt 2 V m t 



Mi 



pi 



-1/2 

< 0.625 (5.24) 



where M Pi ~ 1.2 x 10 19 GeV is the Planck mass. Note that [3\ s has a strong depen- 
dence on A, hence as A(/xo) becomes close to upper bound, the allowed band of Xs(fio) 
for a perturbative positive As(/i) becomes narrower and eventually disappears. 

5.3 Numerical analysis 

The aim of our numerical analysis is to see if the demand of vacuum stability con- 
strains our model parameters which should satisfy EWPT, DM relic density and 
DM direction detection bound, and if there is any way to probe the model in future 
collider experiments. In this regard, the crucial parameters are 

mi, m 2 , a, A, X H , X H s, A s . (5.25) 

These parameters are involved in the following constraints. 

• EWPT: mi, m 2 and a 

• DM relic density: The first 4 parameters 

• DM direct searches: The first 4 parameters 

• Vacuum stability (A#(//) > and As(/i) > 0): All of them. 

Note that not all of those parameters are independent. For example, if we choose mi, 
m 2 and a as an input, each element of the mass matrix is fixed through Eq. (2.7). 
Since % = 246 GeV, A^ is fixed by the first line of Eq. (2.8), but we are free to choose 
Xhs an d Xs- Once (Xhs, -Vs) is chosen, ((jlhs, A*s) i s given in terms of vs by the second 
and third lines of Eq. (2.8), respectively, where we assume fis = 0. Therefore, the free 
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parameters we can use in analyzing the scale dependence of dimensionless couplings 
are 

mi, m 2 , a, A, X H s, X s . (5.26) 

Inspired by the recent results of LHC experiments, we take ni\ = 125 GeV and 
vary vri2 and a in the ranges, 

150 GeV < m 2 < 2 TeV , < a < vr/4. (5.27) 

Since we are interested in parameter space where couplings do not blow up, we 
consider A ranging 

0.01 < A < 0.6. (5.28) 

We find that Xhs blows up if \Hs( m t) > 0.4 even for A^(m t ) = X H M (rn t ) and 
Xs{frit) = X{m t ) = 0. As blows up if Xs(m t ) > 0.26 for X{m t ) = 0.6 and XHs{ m t) = 0. 
If Xhs starts from a negative value, As blows up if XHs( m t) < —0.9 for X{m t ) < 0.6 
and Xs{m t ) < 0.26. Based on this observation, we scan the following ranges of our 
parameters. 

- 0.9 < X HS (m t ) < 0.4, 0.01 < X s (m t ) < 0.26. (5.29) 

Fig. 9 shows a distribution of stable(red dots)/unstable(blue dots) vacua in 
(m2, a) plane for 10 4 randomly chosen parameter sets. Left and right panels are for 
positive and negative Xhs, respectively. In both panels, the upper bound of red/blue 
dots are from the perturbativity constraints which we chose 

A < I- (5.30) 

In case of Xhs > 0(left panel), red dots covers whole region below the perturbativity 
bound and go far below the SM bound. They appear even in the region where mixing 
is quite small. This is because non-zero Xhs can provide a large enough loop effect 
on the running of A# to remove the vacuum instability even if it is upper-bounded to 
a rather small value to avoid blow-up. On the other hand, for Xhs < 0(right panel), 
quite small number of red dots(stable vacua) appear. This is because the constraint 
from unbounded-from-below at large field region removes big chunk of parameter 
space, which is clear from Fig. 10 and 11 where distributions of stable vacua shown 
in Fig. 9 are depicted in (A#s ; As) and (A, As) planes, respectively. 

As an example parameter set for stable vacuum, in Fig. 12 and 13, the RG- 
running of dimensionless couplings are shown for loop-effect only(left) and mixing 
only (right), respectively. Comparing both figures, we observe that, because of in- 
creasing magnitude of Xhs, the loop effect raises up the running of A# at high renor- 
malization scale. Contrary to this, tree-level mixing effect is nearly like a simple 
shifting up of A# that can be also seen clearly in Fig. 14 

The triviality and vacuum stability bound of Higgs mass as a function of the 
renormalization scale is shown in Fig. 14. Dashed lines in the figure corresponds to 
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Figure 9. Vacuum stability in (m2,a)-plane. Left: Xhs > 0. Right: Xhs < 0. Blue 
line is the EWPT bound. Green solid line is the bound for stable vacuum when RG-run 
is SM-like. Green dashed lines correspond to Xh = 0.17,0.20,0.25 from bottom to top. 
Red/blue dots indicate stable/unstable vacuum. 




Figure 10. The distribution of stable vacua appeared in Fig. 9 in (A5,Aj^s , )-plane. 



standard model bounds. Solid lines are bounds in our model. Note that in the left 
panel where only loop effect is included the solid-blue instability line is terminated at 
an intermediate scale of \i. It is because the quartic self-interaction of Higgs decreases 
at low scale but increases at high scale, as can be seen in the right panel of Fig. 12 6 . 

6 As can be seen from Fig. 12, if Xhs is small, A^r(/x) can cross-down zero-point at a scale, but 
cross-up at a higher scale. Instability, however, may not exist if SM-vacuum is the global minimum. 
Even if SM-vacuum is metastable, it is okay as long as the tunneling time is longer than the age of 
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Figure 11. The distribution of stable vacua appeared in Fig. 9 in (A,s,A)-plane. Left: 
Xhs > 0. Right: X HS < 0. 
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Figure 12. RG-running of couplings as a function of renormalization scale for mi = 
125 GeV, 7712 = 500 GeV and a = 0, i.e., no-mixing. Red/blue/green/dashed-blue line 
corresponds to Xh/Xhs/^/^-S- 



As shown in the figure, both of loop effect of Xhs an d tree-level mixing can remove 
the vacuum instability easily in our model though they work differently. 

our universe. Cosmological danger of falling down to wrong vacuum may be also avoidable as long 
as Higgs field were in thermal touch with high temperature radiation background after inflation. 
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Figure 13. RG-running of couplings as a function of renormalization scale for mi = 
125 GeV, iri2 = 500 GeV and a = 0.1, but Xhs = 0, i.e, mixing but no-loop correction. 
Red/blue/green/dashed-blue line corresponds to Xh/ Xhs I X/ As- 




Log[/j/GeV] Log[/i/GcVJ 

Figure 14. The mass bound of SM-like Higgs (mi) as a function of energy scale for 
(a, Xhs) = (0, 0.2)(left),(0.1, 0) (right) with X s = 0.1 and A = 0.4. The red/blue line 
corresponds to triviality /vacuum-stability bound in SM(dashed) and our model(solid). The 
dashed black line corresponds to m-i = 125 GeV. 



5.4 Brief Summary 

In brief summary, the numerical analysis shows that the vacuum stability of Higgs 
potential and perturbativity of couplings constrains new dimensionless couplings of 
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our model as follows. 

< |A| < 0.6 , -0.2 < X HS < 0.4 , < X s < 0.2. (5.31) 

On the other hand, we cannot extract any useful bounds on the 2nd Higgs mass rri2, 
since it depends on an unknown quantity, the singlet VEV Vs- 

6 Conclusion 

In this paper, we have made a careful and comprehensive study of vacuum structure 
and vacuum stability in the singlet fermion dark matter with a real singlet scalar mes- 
senger. We found that the vacuum structure of this model has a very rich structure. 
In the parameter space we explore here, the EW, I, III and SYM phases (definitions 
are given in the text) appear, and any of them can become the global minimum. It 
is remarkable that the region where the EW vacuum is the global minimum is overly 
limited, eliminating the large parameter space. The tree-level analysis shows that 
the global vacuum condition yields the constraint on the mass of the second Higgs 
boson, i.e., \/Xs/2\vs\ < va% < y/2\g\vs\ for v$ ^> vh and a < 0.2 rad. 

At the one-loop level, on top of the SM particles and singlet Higgs, the DM 
also participates in the effective potential. We found that effects of the DM on 
the vacuum structure can be significant depending on the magnitude and sign of A. 
However, it tends to be small once we take = 1112/2 which is favored by the DM 
relic abundance requirement. In such a case, above tree-level bound on ni2 is still 
valid at the one-loop level as verified numerically. Our findings also applies to the 
SM extension with a singlet scalar boson without CDM. 

We also studied vacuum instability caused by the RG-running effect of the top 
quark at high-energy scales. Unlike the SM, Xh in the current model can be positive 
all the way to the Planck scale. The reason is two-fold: due to the tree and loop 
effects coming from the singlet Higgs sector. 

The former is related to the mixing between the doublet and singlet Higgs bosons 
at the tree level. Because of this, the Higgs quartic coupling and mass is no longer 
one-to-one correspondence. We can take larger Xh with mi = 125 GeV by increasing 
a and/or m>2 as shown in Eq. (5.6). This increment in the initial value of Xh can 
avoid zero-crossing up to the Planck scale. 

The latter is nothing but the effect of Xhs on the RG-running of Xh- Such an 
effect can be sizable enough to compensate the negative contribution coming from 
the top quark at high-energy scale, preventing the Higgs potential from generating 
an unbounded-from-below direction or a potentially new global minimum. The per- 
turbativity of the quartic couplings up to the Planck scale are also investigated. The 
results are summarized in Eq. (5.31). 

Before closing, we make a comment on the bound of rri2. Although we obtain the 
strong bound on rri2 through the global vacuum condition, rri2 is still not predictable 
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since Vg and \$ are totally unknown. Therefore, we may need additional information 
to constrain those parameters somehow. For example, since strong first-order elec- 
troweak phase transition as needed for successful electroweak baryogenesis is closely 
related to specific vacuum structure and so the specific vs and As, we may get more 
useful bound on m 2 . We leave this possibility to a future study. 

7 Note Added 

While we are finalizing this work, CMS and ATLAS Collaborations reported a new 
boson of mass around 125 GeV, which might be consistent with the SM Higgs boson. 
The invisible branching ratio of the observed particle seems to be small. In Fig. 15, 
we show the reduction factor for one of the Higgs boson mass is equal to 125 GeV. 
(Note that we took m h = 120 GeV in Ref. [4].) 

Implications of this new results on the model studied in this work are the fol- 
lowing: 

• Since the observed properties of a new boson is close to those of the SM Higgs 
boson, its singlet component should be small, namely the mixing angle a should 
be small in our model. 

• If mi = 125 GeV, the singlet fermion DM mass may have to be greater than 
mi/ 2 ~ 63 GeV so that Hi — > ipip is kinematically forbidden. The current 
Higgs search results can not be directly applied to if 2, since if 2 would be 
mostly a singlet and difficult to be produced at colliders. Also it may decay 
into a pair of DM's with a substantial branching ratio. Note that < 0.3 from 
Fig. 10 of Ref. [4]. 

• If m 2 = 125 GeV, the lighter Higgs boson Hi can be light and may have escaped 
the detection if it is mostly SM singlet. This is still consistent with all the data 
available as of now. 

• In order to test the idea of Higgs portal dark matter with a singlet fermion 
DM, it is important to search for two Higgs-like scalar bosons, one of which 
is to be identified with a new particle with mass around 125 GeV. The other 
scalar could be heavy and have escaped the Higgs search at the LHC and at 
the Tevatron, if it has substantial invisible branching ratio into a pair of DM's. 

• n defined in Eq. (3.4) is always less than one in our model. Therefore, our 
model would be excluded, if turns out to be larger than one (r, > 1) in any 
of the decay channels in the future. On the other hand, if r q < 1 is observed in 
all observed channels in the future, our model could be a good candidate for 
the reason behind it. We have to wait for more precise determinations of for 
all possible measurable decay channels of Higgs-like boson with 125 GeV mass. 
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• In our model, there are two Higgs-like scalar bosons, one of which could be 
mostly singlet scalar and hard to discover at colliders, since r,; < 1. In Fig. 16, 
we show a scattered plot in the (m 2 (i), a)-plane for m 2 > mi = 125 GeV (mi < 
m 2 = 125 GeV). Note that there is no strong constraint on the allowed mass 
range for the 2nd Higgs from the EWPT, collider or DM phenomenology. And 
the current SM Higgs search bounds do not apply directly, since the signal 
strength is reduced in our model (r, < 1) compared with the SM Higgs boson. 
Note that a < 0.4 x 7r/2 for all values of m 2 > 220 GeV, for which case 
H 2 is always singlet-like. We find that r 2 < 0.4 for m 2 > 220 GeV, making 
it difficult to search for the heavier (singlet-like) scalar particle. Therefore 
the current Higgs search should be continued for wider ranges of Higgs mass 
considering a possibility of < 1, in order to look whether another Higgs-like 
scalar boson (mostly singlet-like) exists or not. 

While we are finalizing this paper, there appeared a paper which considers Higgs 
phenomenology of a similar model, a singlet extension of the SM with DM in a hidden 
sector [40]. 
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A The top-quark running mass 

The top-quark pole- mass is related to the MS mass m t (fi) to NNLO as 



M t = m t (fi) 



a 



(n /= 5) 



7T 



(n/=5) 



77 



(A.l) 



where a s nf ^ is the strong coupling with five active flavors, diS are given as [41] 
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with ^ 2 m 1.645 and £3 ~ 1.202 being zeta-constants. For matching, we have to set 
/j, = m t . At this scale the strong coupling can be found from the following relation 
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for the running coupling 



a s (m t ) = a s 
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B Higgs boson masses at the one-loop level 

Here, we write the useful formulae for the Higgs boson masses at the one-loop level. 
Let us begin by listing the field-dependent masses of all particles. The field-dependent 
masses of the Higgs bosons are given by the eigenvalues of the following 2-by-2 mass 
matrix 



M Higgs(^^s) 



— 2 — 2 
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The field-dependent masses of the NG bosons, gauge bosons, top/bottom and singlet 
fermionic DM are respectively given by 
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The one-loop tadpole conditions are cast into the form 
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After imposing the one-loop tadpole conditions, the mass matrix elements of the 
Higgs bosons are 
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Figure 15. Scatter plot in (n,^) plane for the scenario SI, S2 and S3 (from above). The 
points represent 4 different cases: (r2cDM^ 2 ) 3<T , &p (big red), (^cdm^ 2 ) 3ct , (big blue), 
(Qcdm'i 2 ) < j (small orange), and (^cdm^ 2 ) < , &p (small green). (See the text for more 
detail). 
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Figure 16. Scattered plot in (m 2 (i),a) plane for m2 > mi = 125 GeV (mi < mi 
125 GeV). The color scheme is the same with Fig. 15. 
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